Effect of electron-phonon scattering on shot noise in nanoscale junctions 
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We investigate the effect of electron-phonon inelastic scattering on shot noise in nanoscale junc- 
tions in the regime of quasi-ballistic transport. We predict that when the local temperature of the 
junction is larger than its lowest vibrational mode energy eV c , the inelastic contribution to shot 
noise (conductance) increases (decreases) with bias as V (VV). The corresponding Fano factor thus 
increases as yV. We also show that the inelastic contribution to the Fano factor saturates with 
increasing thermal current exchanged between the junction and the bulk electrodes to a value which, 
for V >> V c , is independent of bias. A measurement of shot noise may thus provide information 
about the local temperature and heat dissipation in nanoscale conductors. 
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It is an established fact that for systems with dimen- 
sions much longer than the inelastic mean free path X p h 
(e.g. a macroscopic sample) steady-state zero temper- 
ature current fluctuations (shot noise) are suppressed 
by electron-phonon scattering Q. [3- Q Similarly, for 
metallic diffusive wires with length much smaller than 
Xph (and smaller than the electron-electron scattering 
length), the Fano factor (i.e. the ratio between shot 
noise and its Poisson value, 2el, where e is the elec- 
tron charge and I is the current of the system) equals 
1/3 and is not affected by inelastic processes |4|. Sys- 
tems of nanoscale dimensions may not fall in either 
one of the above cases. In this instance each elec- 
tron, on average, releases only a small fraction of its 
energy to the underlying atomic structure during the 
time it spends in the junction, making transport quasi- 
ballistic H H Q, H H [H EH However, the current 
density and, consequently, the power per atom are much 
larger in the junction compared to the bulk. This leads 
to heating and inelastic features in the differential con- 
duction which are indeed observed in experiments with 
metallic quantum point contacts 113, Iia . Hil Il5| and 
molecular structures || 0, 0, Il7l llSj as a direct con- 
sequence of the interplay between electron and phonon 
statistics |l9|. For these systems it is therefore not ob- 
vious what is the effect of inelastic scattering on shot 
noise. 

In this Letter we show analytically that shot noise in 
quasi-ballistic nanoscale junctions is enhanced by inelas- 
tic scattering whenever electrons have enough energy to 
excite the phonon modes of the junction. The current 
instead decreases. As a consequence, the Fano factor in- 
creases. We find it increases with bias as y/V when the 
local temperature of the junction is larger than its lowest 
vibrational mode energy eV c . We also show that with in- 
creasing thermal current carried away from the junction 
to the bulk electrodes, the inelastic contribution to the 
Fano factor converges to a minimum value independent 
of bias for V » V c . A measurement of the Fano fac- 
tor may thus provide information about the local tem- 
perature and heat dissipation in nanoscale conductors. 



Transport in a model atomic gold point contact will be 
used to illustrate these findings. 

Since the dimensions of the junction are much smaller 
than Xph (and the observed inelastic features in quasi- 
ballistic systems are very small Ell EJ) first-order 
perturbation theory in the electron-phonon coupling cap- 
tures the dominant contribution to inelastic scattering. 
This is the contribution we calculate in this paper. 

Let us assume that the junction is connected to two 
biased bulk electrodes. The electronic states of the full 
system are thus described by the field operator ^ = 



=L,R E 



n (r,K|| 



L( R) 

particle wave functions $ £ 



constructed from the single- 
r, K||) and annihilation 



operators a!^ R ^ corresponding to electrons propagating 
from the left (right) electrode at energ y E . K|| is the 
transverse component of the momentum |21| . We also as- 
sume that the electrons rapidly thcrmalize into the bulk 
electrodes so that their statistics are given by the equi- 
librium Fermi-Dirac distribution, /g = l/(exp[(i? — 
EFL(R))/kBT e ] + 1) in the left (right) electrodes with 
local chemical potential -EVl(.r) , where T e is the elec- 
tronic temperature. In the following we will assume that 
T e = K :20\, and the left electrode is positively bi- 
ased so that Epi, < Efr- The stationary scattering 
states (r, K||) are eigenstates of an effective single- 

particle Hamiltonian H e which may be computed, e.g., 
using a scattering approach within the static density- 
functional theory of many-electron systems [2lJ. The 
combined dynamics of electrons and phonons is described 
by the Hamiltonian (atomic units will be used through- 
out this paper) || 



H = H„+ H, 



ph 



where H ph = | J2 Q 



i,fi£vib 



/ ft 



ph i 



2 2 

to a- 



(1) 

is the phonon 



contribution, with qi^ the normal coordinate and the 
normal frequency of the vibration labeled by the fi-th 
component of the i-th ion. H e ^ p h describes the electron- 
phonon interaction and has the following form 8] 



FIG. 1: Feynman diagrams and corresponding amplitudes 
(see text) of the main electron-phonon scattering mechanisms 
contributing to the correction of the current and noise. 
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where a = L, R and bj V is the phonon annihilation oper- 
ator. {Ain^ju} is the transformation matrix that relates 
Cartesian coordinates to normal coordinates, and J^'^? 
is the electron-phonon coupling constant which can be 
directly calculated from the scattering wave-functions 

J 'eTe 2 = JdrJ dK„*l* (r,K,|) (r,R,)< (r,K„ 

(3) 

where we have chosen to describe the electron-ion in- 
teraction with pseudopotentials V ps (r, R,) for each i-th 
ion [HJ. 

We use as unperturbed states of the full sys- 
tem (electron plus phonon) the states ,Tij U ) 

l^^j (r, K||))® \njv), where rij v is the occupation num- 
ber of the jv-th normal mode. The first-order perturba- 
tion to the wave functions is thus 



\9%*>;n jv ) = \^ W ;n jv ) + \5^ W ;n jv ), (4) 
where the first-order correction term is 



- ujMR). 



\89%;n jv ) = lim E E / dE ' D e' 



a f — L,R j'v' 

(^,;n rv > \H e i- mb \y%; n jv ) ; n rv >) 
e(E,nju) - e(E',rij' V ') - it 



(5) 



with Dg L ' the partial density of states of left (right) 
moving electrons, and e{E,nj U ) = E + [nj V + 1/2) Uj V 
the energy of state \^%; %v) • Carrying out explicitly the 
integrals in Eq. the nonvanishing corrections to the 
wave function can be written as 



FIG. 2: Top panel: ratio of the total conductance G of an 
atomic gold point contact and its value in the absence of 
inelastic effects G° as a function of bias for different val- 
ues of thermal current coefficient (see text): A t h = 10 -19 
(dot), 1CT 17 (dot-dashed), 1CT 15 (dashed), and 00 (solid) 
dyn/(sK 4 ). Bottom panel: corresponding Fano factor ratio. 



+(B? , 2 + Bf vA )\* E _ u . u] n jv - 1), (6) 

where B" vl , B? u2 , B" uZ and B? 4 correspond to the dia- 
grams depicted in Fig. ^ For j^^; rij V }, the coefficients 
are given by: 



r>R 

-°JV,1(2) 



i/j, V ■ >v 

• x /(5 + (n J v))/f(l-/| ±a)i „), (7) 



= and 



r>R 



El 1 A rifJ,,RL 



I / / 



E- u E±uj ji , 



(8) 



where 5 = 1 and "— " sign are for the scattering dia- 
grams (a) and (c); 6 = and "+" sign for diagrams 
(b) and (d). The average number of phonons is given 
by (n,j U ) = 1/ [exp (wjv/fcsTtu) — 11 where T w is the lo- 
cal temperature of the junction [3, E3- Similarly, the 
coefficients in |<5'i , £; nj U ) have the forms 



(9) 



where k = 1, ■ ■ • ,4 ; the notation [L =^ R) means inter- 
change of labels R and L. 
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K the first-order correction to the current is 



i j FR dE j <1R j <1K\ I ',['], 
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(10) 



with lf E = {^ E )*d z {^ E ) - Equation 
ffTTTfl has been simplified by using (i) Ie±u u ,e±w u — 
Ijsei va lid for energies close to the chemical potentials; 
and (ii) I§ B E — —Ie L e , a direct consequence of time- 
reversal symmetry. The current is therefore reduced by 
inelastic effects. 

Let us now calculate the corresponding correction to 
shot noise. We have previously shown that shot noise can 
be written in terms of single-particle scattering states as 

El El 



s 



pEpR 




/ dE 


J dR J dKI^ R E 
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(11) 



which reduces to the well-known formula S oc ^ i Tj(l — 
Ti) when the eigenchannels transmission probabilities Ti 
are extracted from the single-particle states with inde- 
pendent transverse momenta 0, E2, . Replacing 
into ljTT|l we get 



(12) 
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Since the summation over vibrational modes contains 
only positive terms, shot noise is enhanced by electron- 
phonon inelastic effects in the quasi-ballistic regime. 
Therefore, the Fano factor F normalized to the corre- 
sponding value in the absence of electron-phonon inter- 
actions (F°) is 



F/F° = 



(13) 

which increases with electron-phonon scattering. 

Note that due to the orthogonality of phonon states, 
the absolute value of the correction to shot noise is 
smaller than that to the current (cf. Eq. (|10fl and 
Eq. iJUJl). Note also that conservation of energy and the 
Pauli exclusion principle play an important role. The for- 
mer dictates an onset bias V c for inelastic contributions; 
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the latter prohibits the scattering processes depicted in 
Fig. [He) and (d) at T e = K. 

These results are illustrated in Fig. |2 where the inelas- 
tic contribution to the conductance and shot noise are 
plotted for a gold atom placed in the middle of two bulk 
gold electrodes (represented with ideal metals, jcllium 
model, r s w 3). Details of the calculations of both the 
scattering wavefunctions within static density-functional 
theory and the vibrational modes for this system can 
be found in Refs. 0, EI- In the absence of electron- 
phonon interactions, the unperturbed differential con- 
ductance G° is about l-l(m units of 2e 2 /h) and the 
Fano factor is F° ~ 0.14 in the bias range of Fig. H 
Inelastic effects cause a discontinuity in the conductance, 
and a variation of the Fano factor ratio (Eqs. I|13fl ). at a 
bias V c w 11 mV, corresponding to the energy of the 
lowest longitudinal mode of the system. In addition, 
the above inelastic corrections depend on the local tem- 
perature of the junction T w (see Eqs. Q) through (JSJ) 
which, in turn, is the result of the competition between 
the rate of heat generated locally in the nanostructure 
and the thermal current Ith carried away into the bulk 
electrodes H S S E E3 ■ The latter has a tempera- 
ture dependence of Ith = A t hT 4 where the constant 
A t h depends on the details of the coupling between the 
local modes of the junction and the modes of the bulk 
electrodes. At steady state this thermal current has to 
balance the power generated in the nanostructure, which 
is a small fraction of the total power of the circuit (V 
is the bias, R is the resistance) 0>I3- 

The larger A t h, the larger the heat dissipated into the 
bulk and, thus, the lower the local temperature T w |25). 
In the limit of infinite Ath, i.e. T w = 0, at any given bias 
larger than V c , electrons can only emit phonons [{riju) = 
in Eqs. Q and JHJ]. The inelastic contribution to the 
conductance and Fano factor, therefore, saturate to a 
specific value (see Fig. 0). We can derive both the bias 
dependence and this saturation value, to first order in 
the bias, as follows. 

By equating the thermal current Ith to the power gen- 
erated in the junction, it is easy to show that T w — 
ay/V 0) H3 , where the constant a depends on the de- 
tails of the thermal contacts between the junction and 
electrodes. Let us assume for simplicity a single phonon 
mode of frequency lu. For T w > uj/ks, we expand 
{uj V ) w k B T w /to in Eqs. (J7J and ijHJ). From Eq. (|TU|) 
we then get 



G_ 



a- — 7j0(k 



v c )W, 



(14) 



where 8(V — V c ) is the Heaviside function; 7/ = 
\{dI/dV) j (dI°/dV)\ is the relative change in conduc- 
tance due to inelastic effects at V c (its value is about 1% 
for the specific case, in agreement with experiments on 
similar systems 0, 0|). The inelastic contribution to 



4 



the conductance thus decreases with bias as yV. This 
square-root dependence is clear in Fig. |2] for A t h < 10~ 15 
dyn/(sK 4 ) which corresponds to tem pera tures for which 
the condition T w > Lu/kB is satisfied |27j . 

The same analysis applied to shot noise leads to 



S° ~ 



) ~f s 6(V-V c )(V-V c ), (15) 



where 7s = \ (dS/dV) / (dS°/dV) \ is the relative change 
of shot noise due to inelastic effects at V = V c (it is about 
0.04 % for the specific gold quantum point contact). The 
inelastic correction to shot noise thus increases linearly 
with bias for T w > u/ks- Consequently, F/F° cx y/V as 
it is also evident from Fig. [3 

In the opposite limit of perfect heat dissipation in the 
bulk electrodes, i.e. for T w — ► [see Fig. |21 Ath — > 00 
dyn/(sK 4 )], then from Eqs. 10 and © it is easy to prove 
that J/Jq = 1 - 6{V - V c )ji(V - V c )/V and S/S = 
1 + is [{V - V c )/V] 6(V - V c ). Therefore, 



F/F° 



(16) 



which tends to the constant value F/F° — > (1 + 75) /(l — 
7/) as V » V c . 

We have thus shown that the Fano factor depends sen- 
sitively on the efficiency of heat dissipation in nanoscale 
junctions. It therefore provides a tool to probe local tem- 
peratures and heat transport mechanisms in these sys- 
tems. The predictions reported here should be readily 
tested experimentally. 
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